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DIRECT NUMERICAL SIMULATION OF REACTING FLOWS 
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The objectives of thio work are (i) to to extend the technique of direct numeri- 
cal simulations to turbulent, chemically reacting flows, (ii) to test the validity of 

the method by comparing computational results with laboratory data, and (iii) to use 
the simulations to gain a better understanding of the effects of turbulence on chemi- 
cal reactions. In particular, we address the effects of both the large-scale struc- 

ture and the smaller-scale turbulence on the overall reaction rates, examine the rela- 
tionship between infinite reaction rate and finite reaction rate chemistry, and com- 
pare some of the results of our calculations with existing theories and laboratory 
data. The direct numerical simulation method involves the numerical solution of the 
detailed evolution of the complex turbulent velocity and concentration fields. Using 
very efficient numerical methods (e.g., pseudospectral methods), the fully nonlinear 
(possibly low-pass filtered) equations of motion are solved and no closure assumptions 
or turbulence models are used. Statistical data are obtained by performing spatial, 
temporal, and/or ensemble averages over the computed flow fields. 

The scope of work to do this involved the following. First, existing computer 
codes were modified to treat the present problem. Next, extensive numerical testing 
of the computer codes was performed. Finally, in order to examine the effects of the 
mixing layer turbulence, both the large-scale structure and the smaller-scale turbu- 
lence, on the overall reaction rates, a sequence of three problems was computed: (i) 

reactions on a unidirectional (one-dimensional) mixing layer, (ii) reactions on & 
mixing layer experiencing large-scale, two-dimensional vortex rollup, and (iii) reac- 
tions on a three-dimensional turbulent mixing layer. The simulations of the two- 
dimensional mixing layer with vortex rollup are intended to model the large-scale 
structure in the mixing layer, whereas the three-dimensional simulations contain both 
the large-scale structure and the smaller-scale turbulence. 

The numerical testing involved the comparison of computed results with exact solu- 
tions for a number of different cases. Both rigid body rotation and vortex rollup 
flow fields were used. The work greatly extended the results of Orszag (tef. 1) for 
the advection of a passive scalar on a rigidly rotating flow field (the color problem) 
to include also diffusion, chemical reaction, and more complex flows. We have found 
that high accuracy of the spectral methods observed in the advection case is also 
obtained when these further complications are present. Our results indicate that 
spectral numerical methods may prove to be useful in the future both for solving the 
model equations for combustor processes as well as for future studies of chemically 
reacting turbulent flows employing direct numerical simulations* 

The approach of direct numerical simulations allowed extensive examination and 
interpretation of the reaction process. From the one-dimensional simulations , we 
found that we could easily compute finite reaction rates near the fast reaction limit, 
that the results were fairly insensitive to the initial conditions (for the class of 
initial conditions computed), and that the results were in reasonable agreement with 
theoretical predictions. 
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For the two-dimensional simulations, we found that in all of the cases computed, 
the vorticity field and the product field approximately coincided. From the computa- 
tion of volume averages, we observed the enhancement of the overall reaction rate due 
to the vortex rollup. It also appeared that the merging of vortex cores was a more 
significant mechanism in increasing the overall reaction rate than the straining of 
the reaction interface. Significant species segregation was apparent, so that, for 
example, the product of the average concentrations was approximately equal to and 
opposite the correlation between the fluctuating concentrations. Some of the hig.er 
order correlations were also examined. 

In the three-dimensional simulations, the contour plots indicated that the vortex 
rollup takes longer to develop, and the vortices and braids are not as distinct as in 
the two-dimensional case. Also, the vortices that develop are not strongly corre- 
lated laterally. Both the contour plots and the statistical results indicated that 
the spatial segregation was also not as strong as in the two-dimensional case, 
probably due to the weaker vortex rollup as well as the effects of smaller-scale , 
three-d imensional turbulence. 

Comparisons were made between the simulation results and results using similarity 
theory. Approximately linear growth rates of various computed length scales, in- 
cluding the mean velocity half-width, the mean vorticity thickness, and the mean pro- 
duct thickness, were obtained and were in agreement with the theory. Similarity 
scaling was found to collapse quite well the results for the average reactant concen- 
trations, the rms fluctuating reactant concentrations, the concentration correlations, 
the average product concentrations, and the rms fluctuating product concentrations. 

Some limited comparisons were made with laboratory data. Computed profiles that 
were qualitatively similar to corresponding laboratory profiles were obtained for the 
average reactant concentrations, the rms fluctuating reactant concentrations, the 
average product concentrations, the rms fluctuating product concentrations, and the 
concentration correlations. 

We have made some comparisons with existing theories. Donaldson and Hilst have 
suggested, in addition to using the equations for the average concentrations, in- 
cluding the equations for the concentration fluctuations and correlation (ref. 2). 
These equations can be closed by neglecting certain triple moments when compared to 
certain lower order terms. However, our results indicate that the triple moment 
terms are as important as other terms in the equations, so that an assumption of this 
type will probably lead to poor predictions. Mason and Spaulding have proposed a 
model for the mean reaction term, suggesting that it will be proportional to the 
average concentration of the lean species divided by a turbulent time scale (ref. 3). 
We have found that such an assumption will only be moderately successful if applied 
to our case. Finally, Toor has suggested estimating the concentration correlation of 
the reacting species in terms of that for the nonreacting case, which is much easier 
to model (ref. 4). Although this was proposed mainly for statistically homogeneous 
flows, we find that it is a reasonable approximation for our reacting flow 
s imulations . 
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OBJECTIVES 


• Begin the Extension of Direct Numerical Simulations to Combustor 
Flows 

• Develop computational methodology 

• Use Direct Numerical Simulations in Order to 

• Better understand the effects of turbulence on chemical reactions 
. Suggest model improvements 

• Test Validity of the Method by Comparing Computational Results 
with Laboratory Data and Theoretical Models 


PROBLEM CONSIDERED 

• Temporally Growing, Three-Dimensional Mixing Layer 

• Binary, Single-Step, Irreversible Chemical Reaction 

(A + nB -*> products) with negligible heat release 

• Very Small Mach Number 

• Initial Conditions - NonPremlxed Species 







SCOPE OF WORK 


ORIGINAL PAGE fg 
OF POOR QUALITY 



• 



• Modify Existing Computer Codes to T reat the Present Problem 

• Perform Numerical Testing of the Computer Codes 

• Compute the Following Sequence of Problems 

• Unidirectional (one-dimensional) flow 

• Large-scale, two-dimensional vortex rollup 

• Three-dimensional turbulent mixing layer 
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EQUATIONS OF MOTION 


• Conservation of Momentum 
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• Conservation of Mass 
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• Conservation of Species A and B 
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• Initial Boundary Conditions 
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DIRECT NUMERICAL SIMULATIONS 


• Compute the Solution of the Navier-Stokes Equations Numerically 

• Initial Conditions (or Upstream Conditions) 

Random, turbulent-like - so that the flow field starts out 
as turbulent-like - or 

Small amplitude perturbations - compute the flow through 
transition to turbulence 

• Compute the Time Development of the Detailed Structure of the 
Flow Field 

• No Reynolds Averaging 

• Statistical Results are Obtained by Performing Spatial, Temporal, 
or Ensemble Averages over the Computed Flow Field 

• Analogous to Laboratory Experiments 
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NUMERICAL METHODS 


• Pseudo-Spectral Methods 

• Expand the dependent variables in fast-converging series 

E - v 

uM - ikj<K G<k,e ' ‘ 

• Resulting Equations 


JjC^k.t) + i£ u.(k-k')k! fyk') = -R^ C^k-kX^k'l-Dik^lk) 


Reaction 


Diffusion 
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NUMERICAL TESTS 




• Color Problem (Rigid Body Rotation of Passive Scalar) 

• Convection Alone (Orszag, Seinfeld et al.) 

• Convection Plus Reaction (Exact Solutions Available) 

• Convection Plus Diffusion (Exact Solutions Available) 

• One-Dimensional, Diffusion-Reaction Problem 

• Test Against Exact Solutions 

• Two-Dimensional Vortex Rollup 

• Convection, Diffusion, and Reaction 

• "Exact" Solution for 6 = C - C is Known 

from Previous Calculation ^ 2 


NUMERICAL TEST PROBLEMS 


t Two Dimensional Rigid Rotation (Color Problem) 

• Initial condition: C(x,z,0) = exp 

\ *o/ 

• Analytical Solution for Identical Initial Concentration Fields 
(No Diffusion) 

C(x’,z',t) = 

1 + RtC(x',z',0) 

• Analytical Solution with Diffusion (No Reaction) 


C(x',z',t) = 


jFfrf) 


exp 


x' 2 lxl 
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Initial Concentration Field for Color Problea Teats 



Computed Concentration Field After One Revolution - No Diffuaion 
or ReacLion 
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TWO-DIMENSIONAL SIMULATIONS 


• Case I: Most Unstable Mode Based Upon Linear Theory 

(Fundamental) 

• Results in a Single Rollup 

• Case II: Subharmonic of the Most Unstable Mode 

• Results in a Single Rollup, of Approximately Twice the Size 
and Time Period 

• Case III: Fundamental and Subharmonic 90° Out-of-Phase 

• Results in a Double Rollup 
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EQUATIONS (CONT.) ! 

t 

I 

• Define Conserved Scalar 0 = C A - Cg * 

Then 0 satisfies 

Q] -2-A+ u-V0 = D^- Q (No R dependence) 

0t 3x2 j 

with initial conditions 0(t = 0) = 0(0) = C A (0) - Cq(0) ; 

i 

• Infinite Reaction Rate Limit i 

i 

(0 0 *> r _( 0 0*0 | 

1 (o 0<O 2 ( -0 0<O , 

Cp can be obtained from conserved scalar 
0 — Cp + C A 

where 0 also satisfies Equation Q] ; 

with initial conditions 0(t = 0) = 0(0) = C A (0) • 
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CONCLUSIONS 


• Color Problem Tests Showed That 

• Spectral methods produced high accuracy solutions in the 
presence of advection, reaction, and diffusion 

• Weak spatial filter insured numerical stability 

• One Dimensional Simulations Showed That 

• Calculations close to the infinite reaction rate can be done 

• Results were fairly insensitive to initial conditions 

• Results were in reasonable agreement with theoretical 
predictions 

• Two Dimensional Simulations Showed That 

• Reaction zone was significantly stretched by vortex rollup 

• Vorticity and product fields coincided 

• Significant species segregation and "flame shortening" 
effects were present 


CONCLUSIONS (CONT.) 

• Three Dimensional Simulations Showed That 

• Good agreement with self-similarity theory was obtained 

• Technique can treat entire range of reaction rates - from very 
slow to infinite 

• Computed profiles qualitatively similar to profiles from 
laboratory data 

• Quantitative comparisons will be made soon 

• Nondimensional product thickness was 0,22 - compared to 
values of 0.3 to 0.35 from laboratory data 

• No adjustable parameters are used 

• Some presently used estimates of correlation products 
can be significantly in error 


221 



